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Abstract 


Purpose — The purpose of this paper is to investigate the effect of reinforcement, inhomogeneity and initial 
stress on the propagation of shear waves. The problem consists of magneto poroelastic medium sandwiched 
between self-reinforced medium and poroelastic half space. Using Biot’s theory of wave propagation, the 
frequency equation is obtained. 

Design/methodology/approach — Shear wave propagation in magneto poroelastic medium embedded 
between a self-reinforced medium and poroelastic half space is investigated. This particular setup is quite 
possible in the Earth crust. All the three media are assumed to be inhomogeneous under initial stress. The 
significant effects of initial stress and inhomogeneity parameters of individual media have been studied. 
Findings — Phase velocity is computed against wavenumber for various values of self-reinforcement, 
heterogeneity parameter and initial stress. Classical elasticity results are deduced as a particular case of the 
present study. Also in the absence of inhomogeneity and initial stress, frequency equation is discussed. 
Graphical representation is made to exhibit the results. 

Originality/value — Shear wave propagation in magneto poroelastic medium embedded between a self- 
reinforced medium, and poroelastic half space are investigated in presence of initial stress, and inhomogeneity 
parameter. For heterogeneous poroelastic half space, the Whittaker’s solution is obtained. From the numerical 
results, it is observed that heterogeneity parameter, inhomogeneity parameter and reinforcement parameter 
have significant influences on the wave characteristics. In addition, frequency equation is discussed in 
absence of inhomogeneity and initial stress. For the validation purpose, numerical results are also computed 
for a particular case. 


Keywords Shear wave, Self-reinforcement, Magneto poroelastic, Half-space, Inhomogeneity, 
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Paper type Research paper 


1. Introduction 
Study of mechanical behavior of self-reinforced material has a huge importance in the domains 
such as structural engineering and geomechanics. Self-reinforced materials attract many 
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investigators in the field of engineering applications due to their superiority over other 
structural materials in the applications where high strength and stiffness and light weight 
components are required. Most of the structures on or near the surface of the Earth are self- 
reinforced and are constructed in such a manner that the effect of earthquake or similar 
disturbances causes minimum destruction. Many elastic fiber-reinforced composite materials 
are strongly anisotropic in nature. There are many materials for example, the Earth’s crust and 
mantle which are inhomogeneous. The Earth crust contains some hard and soft rocks that 
exhibit both self-reinforcement property and inhomogeneity. As a result of stress developed 
within the crust, deformation takes place along with fracture. These fractures release large 
amount of energy which give rise to elastic waves travel and finally reach to surface. Moreover, 
the Earth crust made of diversity of igneous, metamorphic and sedimentary rocks. These rocks 
are capable to generate magnetic field due to presence of some minerals like iron, nickel and 
cobalt in them. For the said reasons, these rocks can be treated as magneto poroelastic solids. 
An increasing attention is being devoted to the interactions between magnetic and strain fields. 
Belfield et al. (1983) introduced the continuous reinforcement at every point of an elastic solid. 
Vardoulakis (1984) studied the propagation of torsional surface waves in heterogeneous elastic 
media. Propagation of shear waves and transverse surface waves in self—reinforced magneto 
elastic bodies is investigated by Verma (1986) and Verma et al. (1988). Chattopadhyay and 
Chaudhary (1995) discussed the propagation of magneto elastic shear waves in an infinite self- 
reinforced plate. Chaudhary et al. (2006) reveal that reflection and transmission coefficients are 
strongly influenced by the reinforcement parameters, and angle of incidence of the incident 
wave. Chattopadhyay ef al (2010) investigated the dispersion of shear wave in multilayered 
magneto elastic self-reinforced media. Effects of point source, self-reinforcement and 
heterogeneity on the propagation of magneto elastic shear waves are presented in 
Chattopadhyay et al (2011). Samal and Chattaraj (2011) discussed the surface wave 
propagation in fiber-reinforced anisotropic elastic layer between liquid saturated porous half 
space and uniform liquid layer. Propagation of shear waves in an irregular self-reinforcement 
layer, the combined effects of reinforcement, magnetic field and irregularity are studied by 
Chattopadhyay and Kumar (Chattopadhyay and Singh, 2012). Chattopadhyay et al (2012) 
discussed the torsional surface waves in a self-reinforced medium over a heterogeneous half 
space. Using the dual reciprocity boundary element method (DRBEM), Fahmy (2012a, 2012b, 
2012c, 2012d) studied the transient magneto-thermo-viscoelastic stresses in a non- 
homogeneous anisotropic solid subjected to a heat source. The effect of a moving heat source 
on the transient magneto-thermo-visco-elastic stresses in a non-homogeneous anisotropic solid 
is investigated (Fahmy, 2012a, 2012b, 2012c, 2012d). Fahmy (2012a, 2012b, 2012c, 2012d) 
discussed the generalized magneto thermoelasticity problem of rotating anisotropic viscoelastic 
functionally graded solids. Generalized magneto-thermo viscoelastic problem of rotating 
functionally graded materials are investigated by Fahmy (2012a, 2012b, 2012c, 2012d) and 
Fahmy (2013). Abd-Alla ef al studied the surface waves propagation in fiber-reinforced 
anisotropic elastic media subjected to gravitational field (Abd-Alla et al, 2013). Chattaraj and 
Samal studied the propagation of Love waves in fiber-reinforced layer lying over a gravitating 
anisotropic porous half space (Chattaraj and Samal, 2013). From Chattaraj and Samal (2013), 
authors concluded that the transverse and longitudinal rigidity of reinforced material, and 
porosity of half space have a significant effect on wave propagation. Chatterjee ef al. (2015) 
studied the propagation of shear waves in visco-elastic heterogeneous layer overlying on 
initially stressed half space. Effect of reinforcement on half-space is discussed (Singh et al, 
2016). Tong et al. (2016) investigated wave propagation characteristics in fluid saturated porous 
materials in the framework of nonlocal Biot’s theory. By considering the interface between self- 
reinforced medium and half space, effect of reinforcement and heterogeneity on the propagation 


of torsional surface waves is discussed by Paswan et al. (2017). Effect of gravity and initial 
stress on propagation of torsional surface waves in heterogeneous medium is studied 
(Mukhopadhyay ef al, 2017). The reflection of plane waves in a rotating magneto elastic fibre- 
reinforced semi space is investigated by Roy et al. (2017). From Roy et al. (2017), it is clear that 
magnetic field, rotation, fibre-reinforcement and surface stress play vital role in wave 
propagation, particularly, in reflection of plane waves. Using Biot’s theory (Biot, 1956), shear 
wave propagation in magneto poroelastic dissipative isotropic medium sandwiched between 
two poroelastic half spaces is studied (Balu et al, 2017). In Balu et al (2017), wave 
characteristics are computed at the fixed heterogeneous and magnetic poroelastic coupling 
factor. Tong et al. (2018) studied the Rayleigh wave propagation in porous fluid saturated 
materials and concluded that the nonlocal parameter does not have significant influence on the 
characteristics of Rayleigh waves within a low frequency range when compared to that of 
classical Biot theory. Effect of several parameters of reinforced sandwich porous plates such as 
aspect ratios, volume fraction, types of reinforcement and thickness of plate on bending is 
investigated by Medani et al. (2019). The nonlinear static, buckling and vibration analysis of 
viscoelastic micro-composite beam reinforced by various distributions of boron nitride 
nanotube (BNNT) is studied by using finite element method (FEM) (Alimirzaei et al., 2019). 
Vibration analysis of functionally graded micro beams is studied by Tlidji et al. (2019). Boukhlif 
et al. (2019) made dynamic investigation of functionally graded (FG) plates resting on an elastic 
foundation using a simple quasi-3D HSPT in which the stretching effect is considered. Some 
recent papers (Karami et al., 2019a, 2019b; Karami et al., 2019a, 2019b) are reported in the field 
of the functionally graded nanoplates reinforced with graphene nanoplates and nanoshell in 
magnetic field. In the above-mentioned papers, the authors have not studied the wave 
propagation in magneto poroelastic medium sandwiched between self-reinforced medium and 
poroelastic half space. Hence, in the present work, the same is investigated. This particular 
setup is quite possible in the Earth’s crust. All the three media are assumed to be 
inhomogeneous under initial stress. The significant effects of initial stress, and inhomogeneity 
parameters of the individual media have been studied. 

This paper is organized as follows. In Section 2, formulation of the problem is given. The 
solution of the problem is discussed in Section 3. Boundary conditions and frequency 
equation are presented in Section 4. In Section 5, particular cases are discussed. Numerical 
results are described in Section 6. Finally, conclusion is given in Section 7. 


2. Formulation of the problem 

Consider the magneto poroelastic medium sandwiched between inhomogeneous self-reinforced 
medium and poroelastic half space, all assumed to be under initial stress. The origin O is 
considered in a manner such that it separates the sandwiched layer from upper and lower half- 
spaces. The horizontal axis is assumed to be along 7-axis, and z-axis is in downward direction 
as shown in Figure 1. From the figure, it is clear that upper self-reinforcement layer M4 (say), 
magneto poroelastic layer M (say), and poroelastic half space 3 (say) occupy the spaces 
(hy < z < hy), < z < 0) and & > 0), respectively. Equations of motion for shear 
deformation in the isotropic poroelastic medium in absence of dissipation are (Biot, 1956) 
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In equation (1), v and V are the displacement components of solid and fluid, respectively, 


P=A + 2N,Q, and R are all poroelastic constants, e and « are the dilatations of solid and 
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Figure 1. 
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fluid. The mass coefficients p11,p12 and pə satisfy the relations pı = p11 + p12 =U-B) ps, 
p2= p12 + p2 = B pf p = pı + p2=ps + B(pe— ps). Here, pı and p2 are masses of solid 
and fluid per unit volume of aggregate, p is total mass of solid-fluid aggregate per unit 
volume. p, and pyare mass densities of the solid and fluid, respectively.@ is the porosity of 
the aggregate. The solid stresses oj, and fluid pressure s are given by (Biot, 1956): 


oy = 2Ne;j + (Ae + Qe) dy, (2,7 = 1,2,3), s=Qe+Re. (2) 


In equation (2), es are strain components, and ô; is the well-known Kronecker delta 
function. 


3. Solution of the problem 
The solutions for three medium are presented in the following sub sections. 


3.1 Upper self-reinforced medium(-h»y < z < —h;) 

Let (u,0;,W;) and (U;,,V;,W,) be the displacement components of solid and fluid, 
respectively, in the upper self-reinforced medium. For the shear vibrations, displacement 
components are “1 = w1 = 0, vı = 01(7,z,0), and U; = W, = 0, Vi = Vi(7,z,). The stress tensor 
for a self-reinforced medium along the preferred direction @’ is given below (Chattopadhyay 
et al., 2011). 


Tij = Depp d jj +2u TE + a(araer 6 ij + illjekk) + 2( My — UT) (Gidpej + Gj Api ) 
a BApAinCem GG; ’ 
Here, 7/,k,4m = 1,2,3, and the coefficients A,a,8,47,41 are the elastic constants with the 


dimension of stress,u7 and mw, are the transverse and longitudinal shear modulus, 
respectively, the vector 7 = (a1, a2, a3) T is the preferred unit direction of reinforcement so 


that a? + a2 + a2 = 1. In the plane under consideration, 7; = (sin¢,0,cos@)", ¢ is the 
angle of inclination to the radial axis. For torsional vibrations, the strain components are as 
follows: 


0 ue ande = 3 (2) 
Cry 200 Ezz Ezr > 62 2 oz’ rọ — 2 s 


Or r 
Using the above in equation (2), the non-zero stress components are given by: 
w v Ov Ov v Ov 
Oro = LT fs; (2 — 2) + s% 092 = HUT Is = -= *) + s2 ; (3) 


In equation (3), Sı = (1 + (u — ja): So = (u'-1) as, S= 1+ (u — 1)az, 
w= er = BL, UT = URE”. 
In this case, the equations of motion are: 


0 ð 2 0 oO 
r H z H r P 2) = { y; ij 
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o 
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Because of the inhomogeneity with respect to the z coordinate the following are assumed in 
equation (4), P;denotes the initial stress, and P, = Poe”, p11 = P110”, P12 = P120”, P22 = 
P2202”, here a is a constant, P10,P110, P120 P220, are the values of their respective parameters 
atz = 0 For harmonic waves, the displacement component v(z) is assumed as: 


v (2) = Akri, 5) 
Here, k,c,w(=kc) and j;(.) are wavenumber, wave velocity, circular frequency and Bessel’s 


function of first order and first kind, respectively. From equations (3), (4) and (5), the 
following equation is obtained: 
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v1 (2) = exp(—Ez/2) (Cicos( x12) + Cosin( x12) i (kr)?! (7) 
In equation (7), x; =F— ey , C1, C2 are arbitrary constants. The pertinent stress component 
Toz 1S given by: 


Ji (kr) 
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3.2 The sandwiched magneto poroelastic medium(—h; < z < 0) 
Let (u2,V2,W2) and (U>, V2, W2) be the solid and fluid displacement components, respectively, 
in the sandwiched magneto poroelastic medium. For shear wave, equation (1) takes the 
following form: 
OTe i 1 OT 96 i OO 62 : 2076 Ow, i (7 B) g 

5 oo + az tp TPT V XB) 4 = ggu + PiwV2), 


0 
0= JE (P1200%2 + P2200 V2)- 
(8) 


In equation (8), (7 x; B) i is the 6 component of electromagnetic force, T is the electric 
current density and B is the magnetic induction vector. The Maxwell’s equations of the 
electromagnetic field are (Balu et al., 2017): 


as 
V-B=0,V xB =O Ux HAT Ban =0(F+ 4B), 


(9) 


here, E is the induced electric field, H is the magnetic field consisting both primary and 
induced magnetic field, u and ø are the induced permeability and conduction coefficient, 


M; 
respectively. The Maxwell stress tensor (7?) -= pe(Hip; — Hipi — Hyppd yj). Let 


— 


H = (H,, Hə , Hz) and change in magnetic field be (P4,P2,P3) If the displacement current is 
absent, then equation (9) reduces to: 


= 
Axial component of equation (10) 
OH, fi] o ee) 0 ( eee) 
= A ' A ll 
Ot or Ç Ot Oz a Oz aa) 
It is assumed that the primary magnetic field is uniform throughout the space. It is clear that 
there is no perturbation in H, and H,, but there is perturbation ps in Hy Then H, = Ho, Ho = 


Hoz + po, and H; = Hog where Ho1,Ho2, Hoz are components of the initial magnetic field H o- 
=> 
As H = (Aocos¢ , p2, Hsing), where Ho =|Ho| and œ is the angle at which wave 


crosses the magnetic field. Thus, H canbe expressed as: 
H = (Hycos$ , po, Hosing ). 


Here, initially it is assumed that p» is zero. Then equations (11) and (12) give: 
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Integration of above equation with respect to n gives: 


po= Hycosp 24 Hsing = = 


(12) 


(13) 


Equation (9) and the relation ees = AV x A, x H+ (A.V). give the 


electromagnetic force as given below: 
(7 x B) = Me a 


In this case, the components of ( J xB ') are 
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For shear harmonic waves, the displacement component v»(z) can be assumed as: 


v9(2) = RDA lkr e”. 
Equations (8), (14) and (15) give: 


Ëh -p d&dh 
d2 | dz dz 


+ Dəf = 0. 
In equation (16), 
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Here, y> = DŽ, C3,Cyare arbitrary constants. Substitution of above displacement 
component in equation (2) gives the following stress component in the medium under 
consideration: 


(Foz), = Ca (woot at: g2)’exp(—Dyz/2)) (cos( w22)x2 = (Dy/2)sin(x22)h(0)) ) 


+C ( (xe + q2)?exp(-Dz/2)) (sin(x22)x2 — (D1 /2)cos( xah) ) 


3.3 The lower poroelastic medium (z > 0) 

In the lower heterogeneous poroelastic half space, let (u3,v3,w3) and (U3, V3, W3) be the solid 
and fluid displacement components, respectively. For shear wave in this space, the equation 
of motion (1) takes the following form 


avs 10v3 v3 Ə Ov3 fod 
ao zo 2) T G3 a) dE (P1100083 + P12000 V3); 
2 


ð 
0= OE (P 1200003 + P2200 V3)- 


(18) 


In equation (18), w3 = p3o(1 + h2), G3 = uall — P3/2), P3 = P3o(1 + laz), €30 = P30/2 430. 11000 
= puol] + 42), P1200 = P1200011 + 42), p22000 = P2200011 + 42). M30, P30, P1100015 
P120001,P 220001, are the initial values of the respective parameters. The non-zero stress 


components area) = 43(Z) (% — ), O29 = b3(Z) (2) . As in earlier cases, the solution 


of equation (18) v3(z) assumed as: 
03(2) = fla yi(kne@’. (19) 


Equations (18) and (19) give: 
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In equation (20),8; = (a E 


To make nonlinear equation (20) linear, assume that 


Fa0(2) 
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here /30(z) is an arbitrary function, G39 = 30(1 — €30), 4 = l1 — E 303/1 — E30, Then equation 
(20) reduces to: 


dfs 
dz 


(2) = (21) 
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substituting f30(z) = (n), equation (22) becomes: 


R 1 
p" (n)ṣ4 (= | - 7) o(m) =0 (23) 


In equation (23), R = me A ( a 1) | x3 X35 li (h — alk equation (23) is the well- 


known form of Whittaker equation (34). The solution of equation (23) ¢ (n) = CsWro(n) + 
CeW_ro-n). Wro, W_ro are Whittaker functions, Cs,C, are arbitrary constants. For 
bounded solution, Cg must be zero. In this case, the solution of equation (20) is: 


_ 2xX3k3 4 )| 1 iot 
hlz) = o |wea( | (1 + hz) Gal tla 1(kr)e (24) 


Using equations (24) and (19) in equation (2), the pertinent stress is given by: 


(o02)y4, = Cs | (a/Gn( + he)"”) | (4(R — 2/2) /(Gn(1 + ua”) 
WR, 0(2y3h3(1 + uz) /lu) + (R — 1/2)? Wr_10(2x3h3(1 + Lz) /ly)] 


— (t4/2G50(1 + tuz)*””) Wao (2xsha(1 + lz) /ls)] 


4. Boundary conditions and frequency equation 
The boundary conditions at free boundary and two interfaces are prescribed as follows: 
Assume that z = —hs is stress free if (a9), = 0, at z = —h;, the displacement and stress 
must be continuous, Le. (Fez), = (Foz) Myo Yl = Y2 and at the plane of contact of the 
magneto poroelastic medium P) poroelastic half space, i.e. at z = 0 the displacement and 
stress must be continuous, Le. (gz) M: = (762) My? ¥2 = 03. These boundary conditions lead 
to the following system of homogeneous equations: 


[Ai] [G] = [Olz 22 =1,2,3,4,5. (25) 


Equation (25) results in a system of five homogeneous equations in five arbitrary constants 
C,,C2,C3,Cy,Cs. To obtain a non-trivial solution, the determinant of coefficients must be zero. 
Accordingly, the following frequency equation is obtained. 


|A; | = 0,7,j= 1,2, 3; 4, 5; (26) 


5. Particular cases 

Case (i): In the sandwiched medium if the magnetic and porous properties are absent, then 
the problem is reduced to that of classical elasticity (Paswan ef al., 2017), and the frequency 
equation (26) is: 
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Case (ii): When the inhomogeneity parameter is absent, i.e. for upper layer a > 0,47 = ur = 


y, for intermediate layer q — 0, and for lower poroelastic half space 4 —> 0, l — 0, l — 0, 
then the frequency equation (26) reduces to: 


0 Aynsin(xuh2) 0 
0 Aysin( xA) Ag, s (28) 
Azcos( ¥>h1) Azisin( xh) As 


In equation (28), 
Au = xk), 
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Aa = u€ ON? (y2-Dyhy/2)hi (hr), Asi = ilk), An = e” hkr. 

H20D, are similar expressions as equation (17) with q—0. 

Case (ii): Apart from the conditions of case (ii), if the initial stress is absent, i.e. for upper 
layer Pı — 0,€19-0, for intermediate layer P> — 0,€29 — 0,Ga9 H2, and for lower 
poroelastic half space P3 — 0, €39 > 0,G39 > [30,k3 —> k. 

Then the frequency equation (26) reduces to: 


0 Bysin(xi2h2) O 
0 Bysin( x 12/1) Ao =0. (29) 
Ascos( xoh) Assin( xəhı) As 


In equation (29), 
By2= x 2/ilkr), 


ý _ Ji (kr) Jilkr)k 1 La? P1100P2200 — P1200 i 
2O hk) he) P2090 ui (kr) 


The other expressions for equation (29) are similar to that of equation (28) with P.—0. 


6. Numerical results 
For numerical process, the following materials are used. 

The material parameters for upper self-reinforced medium (Chattaraj and Samal, 2013) 
(Say Mat-l), are: 


uT = 2.46x10°N /m?, wy, = 5.66x10°N /m?, p11 = 1.92637x 10°?kg /m? 


P19 = -0.002137x10°kg/m?, po = 0.215337x10°kg/m®* 


The intermediate magneto poroelastic medium is assumed to be sandstone saturated with 
kerosene (Yew and Jogi, 1976) (Say Mat-II). The values in this case are mọ = 70199 x 
10°°N/m?, p 1400 = 24152 x 10°kg/m”, p120 = -300kg/m?, P2209 = 500 x 103kg/m?. 

The lower half-space is assumed to be sandstone saturated with water (Fatt, 1959) (Sa 
Mat-IIl). The values in this case are uso = 0.92 x 10’°N/m”, p11000 = 1.9002 x 10°kg/m”, 
P 12000 = 0, P 22000 = 0.2268 x 10kg/m?. 

Using these values in the frequency equation, the implicit relation between the phase 
velocity and wavenumber is obtained. Phase velocity is computed against wavenumber. 
The magneto-poroelastic coupling factor is taken to be u,H?/N = 0.01 (Balu et al., 2017). 
Moreover, the values for #,kr are taken to be 7r/6,5 (Chattopadhyay ef al., 2012) and the 
other parameter values are in various cases specified in Table 1. The values are computed 
using the bisection method implemented in MATLAB, and the results are depicted in 
Figures 2-6. 


ah hh bh Ishy qh €10 &20 E30 h/h bl i 
Figure 2 0.4 0.5 0.4 0.3 0.4 0.4 0.3 0.2 -= 12 
Figure 3 = 0.5 0.4 0.3 0.4 0.4 0.3 0.2 1.6 1.2 
Figure 4 = 0.5 0.4 0.3 0.4 0.4 0.3 0.2 16 — 
Figure 5 0.4 0.5 0.4 0.3 - 0.3 0.2 16 12 
Figure 6 0.4 0.5 0.4 0.3 0.4 -= 0.3 0.2 16 12 


Phase velocity 
feaa pai ae ae eee 
RSE NWRUADATMOON 


15 2 25 3 35 4 45 5 55 
Wavenumber (kh1) 


Magneto 
poroelastic 
medium 


Table 1. 
Parameter values 


Figure 2. 
Variation of phase 
velocity with 
wavenumber for 
various values (u z/ 
u 7) in upper 
reinforcement 
medium 


EC 


Figure 3. 

Variation of phase 
velocity with 
wavenumber for 
various ratios (h/h) 


Figure 4. 
Variation of phase 
velocity with 
wavenumber for 
various values of 
inhomogeneity 
parameter (ahı) in 
upper reinforcement 
medium 


Figure 5. 
Variation of phase 
velocity with 
wavenumber for 
various values of 
inhomogeneity 
parameter (q/,) in 
magneto poroelastic 
medium 


Figure 6. 
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Figure 2 depicts the variation of phase velocity with wavenumber for different values of 
(uL/ ur) in upper reinforcement medium. From this figure, it is clear that, as wavenumber 
increases, phase velocity, in general, decreases. The notations M1/Mt = 1, M1/Mt = 2, M1/Mt 
= 3 used in the Figure 2 represent the different values of(u}, u^). Figure 3 depicts the 
variation of phase velocity with wavenumber for different values of(/2/h1). From Figure 3, it 
is seen that as wavenumber increases, phase velocity, in general, decreases. Figure 4 depicts 
the variation of phase velocity with wavenumber for various inhomogeneity parameter 
values in upper reinforcement medium. From this figure, it is observed that, as wavenumber 
increases, phase velocity, in general, decreases for all values of inhomogeneity parameter 
(ahı). Figure 5 depicts the variation of phase velocity with wavenumber for various 
inhomogeneity parameter values in magneto poroelastic medium. From this figure, it is clear 
that, as wavenumber increases, phase velocity, in general, increases when the 
inhomogeneity parameter gh, = 0.400, whereas when gh, = 0.406,0.412, phase velocity, in 
general, decreases. For the validation, the results are computed for particular case (i). Phase 
velocity is computed against the wavenumber for various heterogeneity parameters in 
presence of reinforcement and free from reinforcement. The results are presented in 
Figure 6. In Figure 6, Curves 1, 2 and 3 pertaining to presence of reinforcement, whereas the 
Curves 4, 5 and 6 pertaining to free from reinforcement for the steel (Paswan et al., 2017). 
From this figure, it is clear that as wavenumber increases phase velocity decreases as the 
values of heterogeneity parameter (10) increases in both reinforced and reinforced free 
cases. These results are in agreement with that of the paper (Paswan et al., 2017). 


7. Conclusion 

Shear wave propagation in magneto poroelastic medium embedded between a self- 
reinforced medium and poroelastic half space are investigated in presence of initial stress, 
and inhomogeneity parameter. For heterogeneous poroelastic half space, the Whittaker’s 
solution is obtained. From the numerical results, it is observed that heterogeneity parameter, 
inhomogeneity parameter and reinforcement parameter have significant influences on the 
wave characteristics. In addition, frequency equation is discussed in absence of 
inhomogeneity and initial stress. For the validation purpose, numerical results are also 
computed for a particular case. 
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